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Abstract 

We characterize the transport properties of functionalized graphene nanoribbons using ex- 
tensive first-principles calculations based on density functional theory (DFT) that encompass 
both monovalent and divalent ligands, hydrogenated defects and vacancies. We find that the 
edge metallic states ai^e preserved under a variety of chemical environments, while bulk con- 
ducting channels can be easily destroyed by either hydrogenation or ion or electron beams, 
resulting in devices that can exhibit spin conductance polarization close to unity. 
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Graphene is currently the subject of intense experimental and theoretical investigations, thanks 
to remarkable and novel physical features .->^>2jii^ Linear dispersion around the Fermi energy 
{Ep), high crystallinity and mobility, ballistic transport on the sub-micrometre scale (even at room 
temperature), and a two-dimensional structure that is amenable to lithographic techniques-^>^ offer 
great promise for nanoelectronics applications, from field effect transistors to interconnects- to 
spintronics and related applications.-^*^ 

Electronic-structure calculations have shown that (real-space) edge states emerge in zigzag 
graphene nanoribbons (ZGNR);— the resulting peak in the density of states (DOS) at Ep induces 
a magnetic instability that leads to a ground state with antiferromagnetic (AFM) ordering, where 
the spins on the two ribbon edges— i^^*^ have opposite orientations. Under an applied electric field 
either half-metallic— »^ or half-semiconducting— behavior can be observed, of extreme interest 
for spintronics applications.— However, the magnitude of the field perpendicular to the ribbon 
required to close the band gap for one of the spins can be as large as 0.1 V A^^ although it is 
a decreasing function of the ribbon width. In contrast with electronic- structure predictions,— re- 
cent measurements'^ show a weak dependence of the energy gap on the crystallographic direction, 
highlighting the critical role that edge structures,— passivations i^iZO ,!! scattering taking place 
at rough boundaries—"^ can play. 

Since edge terminations, passivations and defects play a key role in the performance of graphene 
devices,— we characterize in this Letter the avenues available to control electronic and spin 
transport with chemical functionalizations or ion- or electron-beam treatments. We study single- 
and double-bonded moieties (H, OH, F, CI, Br, S, O, NO2, and NH2) and hydrogenated defects and 
vacancies with extensive DFT first-principles calculations.— We use an original approach suitable 
to describe the electronic structure and ballistic transport in large-scale nano structures ,—>21 based 
on the chemically accurate and minimal basis of maximally-localized Wannier functions (ML- 
WFs).^'29^ 



First, we consider in Figure 1 the band structure and DOS of the ground state for the represen- 
tative single- and double-bonded ribbons H:ZGNR-8, S:ZGNR-8 and 0:ZGNR-8 (we use here the 
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Figure 1 : Band structure and DOS for the AFM ground state of (a) H:ZGNR-8, (b) S:ZGNR-8 and 
(c) 0:ZGNR-8. Dashed and dotted lines and labels identify some relevant band edges (see text). 
The energy scale is such that Ef=0 eV. 

convention of indicating the species functionalizing the edge of the GNR first, and the width of the 
ribbon in graphene units last). As widely reported, the ground state for H:ZGNR-A^ is AFM and 
insulating, with the two bands close to Ef originating from bonding and antibonding Bloch sums 
of the carbon orbitals. A similar description applies to other monovalent passivating species, 
such as OH or F, or even to more complex electropositive or electronegative ligands, such as NH2 
and NO2, since in these latter cases the bands originating from the lone pairs remain below the top 



of the valence. We label the band edges generated by the orbitals A 1 , A2a,b and A3 (see |Figure 1[ ). 
As k goes from the Brillouin-zone (BZ) center to the edge, the states in these two bands evolve 
from being delocalized in real space inside the ribbon (i.e. in the "bulk") to becoming localized at 
the two ribbon edges. As it will be shown later, such observation is critical to engineer the balance 
between unpolarized bulk ribbon conductance and the conductance channels at the ribbon edges. 



Double-bonded ligands, such as O and S (Figure 1 ^b) and (c)), give rise to a different picture: 
the orbitals of the passivating species undergo sp^ hybridization, with one sp^ orbital contribut- 
ing to a a carbon-ligand bond and the other two accommodating two lone pairs; the remaining 
p^ orbital participates in the n carbon-ligand bond. Remarkably, the semiconducting character 
of the hydrogenated ZGNR is reverted to metallic in presence of these edge passivations. In- 
deed, for these double-bonded functionalizations the Ai —A2 and A3 — A2 bands merge even in 
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Figure 2: Band structure (red dots) and its MLWFs interpolation (solid black lines) for 0:ZGNR-8, 
obtained by including, respectively: (a) all MLWFs (see text), (b) all MLWFs but the MLWFs 
on the oxygens, (c) all MLWFs, (d) all o MLWFs, and (e) all o MLWFs and the lone-pair sp^ 
MLFWs. NM state is chosen for clarity. 

the AFM ground state, but at energies higher than Ef\ this can be contrasted with non-magnetic 
(NM) H-ZGNRs, where the bands merge exactly at Ef. The fully occupied A i —Aia band of AFM 
H:ZGNR-8 becomes a mostly occupied band for S passivations ( Figure Ifc )) and a mostly empty 



band for O passivations (Figure 1 c)) - the more electronegative oxygen is extracting charge from 



the ribbon. The double-bonded moieties show two new bands (whose extrema at = are denoted 
by LP[ and LP'^) that remain very close in energy across most of the BZ. Direct inspection of their 
charge density shows a marked lone-pair character, and so they are the analogue of the lone-pair 
bands alluded to before for other monovalent functionalizations (e.g. NH2), with the key differ- 
ence that these now cross the Fermi level. Differences between the AFM ground state and the NM 
solution are much less pronounced for these double-bonded moieties than for H:ZGNR-A^ (see also 
later. ITable II) . the only effect being a small splitting arising for these two lone-pair bands. Indeed, 
the origin of the magnetic instability, arising in H:ZGNR-A/^ from the large DOS of the NM state at 



Ef, disappears in the presence of double-bonded species (Figure 1 'b), (c)). Charge density plots 



for the A1-A2-A3 bands reveal, for 0:ZGNRs, an opposite behavior with respect to hydrogenated 
ribbons: eigenstates near = are mostly localized onto edge O atoms, whereas localization in 
the bulk of the ribbons is found towards the zone boundary. 

Clear chemical insight on the nature of all the energy bands is obtained after we extract ML- 
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WFs— '^'^'^ from the occupied and the n unoccupied manifolds of these ribbons; we choose here 
as an example the NM state of 0:ZGNR-8. The resulting MLWFs map the Bloch bands onto an 
explicit "tight-binding" basis of localized orbitals, composed of one p^^ MLWF for every C and for 
every O, one o bonding MLWF for every C-C and for every C=0 bond, and two lone-pair sp^ 
MLWFs for every O. The band structure of the ribbon in the full BZ can then be straightforwardly 
interpolated across the BZ by diagonalizing the Hamiltonian in this minimal basis set,— >^ or in 



different subsets. We show in Figure 2 our results: the first panel compares the band structure 
obtained by diagonalization in the minimal MLWFs basis and the full calculation in a complete 
plane-wave basis set, showing that the manifolds of interest are all reproduced with extreme ac- 
curacy. The central role played by the oxygen p^ MLWFs is highlighted in panel (b): removing 
these from the minimal basis is sufficient to destroy the agreement around Ef, resulting instead in 
energy bands similar to those of NM H:ZGNR-8, where the Ai — A2 and A3 — A2 bands merge, in 
proximity of the BZ boundary, at Ef. The manifold obtained by all the p^ (including those on the 
oxygens) is shown in panel (c), while the o manifold is in panel (d); in this latter case, quantitative 
agreement is obtained only after the lone-pair sp^ MLWFs on the oxygens are added to the a basis 
set (panel (e)), since these lone pairs not only give rise to the two lone-pair bands LP\ and LP2, 
but are needed to bring the o manifold in full agreement with the exact reference result. The two 
manifolds in (c) and (e) represent a virtually exact decomposition of the original band structure, 
highlighting once again the role of MLWFs in capturing a faithful physical picture in the smallest 
possible representation. 

The two lone-pair bands, originating from the sp^ MLWFs, are very dispersive, with a band- 
width of 3 eV for 0:ZGNR-8 and 9 eV for S:ZGNR-8; the larger size of the sulphur is responsible 
for the broader dispersion, with the largest hopping element in the Hamiltonian being -0.9 eV for 
0:ZGNR-4, and -2.94 eV for S:ZGNR-4. Nevertheless, there is no instability toward dimerization 
at the edges, even if the bond lengths for 0-0 and S-S of 1.48 and 2.05 A are comparable to the 
2.46 A spacing along the ribbon edge. Only halogenation with larger species (CI and Br) gives rise 
to unit-cell doubling, with neighboring halogens assuming opposite tilts with respect to the ribbon 
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Table 1: Unit-cell ground-state formation energy {E^ ), magnetic instability (E^m ~ E^fm), 
magnetic interaction strength (EfM — Eafm) and magnetic moments in the FM state (IIfm) for 
ribbons of different widths and passivations. Energies are in eV, magnetic moments in Bohr mag- 
netons. 



ribbon 




Enm — Eafm 


Efm — Eafm 


IJ-FM 


0H:ZGNR-4 


-3.38 


0.036 


0.011 


0.37 


F:ZGNR-4 


-2.70 


0.047 


0.014 


0.35 


H:ZGNR-4 


0.38 


0.059 


0.014 


0.39 


H:ZGNR-6 


0.42 


0.081 


0.015 


0.45 


H:ZGNR-8 


0.42 


0.083 


0.007 


0.62 


0:ZGNR-8 


-1.03 


0.008 


0.000 


0.19 


S:ZGNR-8 


0.99 


0.002 


0.000 


0.17 



plane. 

The thermodynamic stability of the different ribbons is summarized in ITable 1[ negative for- 
mation energies are found for OH, F and O passivations. The AFM state is the ground state 
for all the systems considered, but the energy difference between the FM and AFM states de- 
creases with the ribbon width (see ITable II and Ref.— ). Room-temperature FM coupling be- 
tween the edges could be stabilized by an applied magnetic field or induced by defects or ad- 
sorbates;— adsorbates can also break the spin-up/spin-down symmetry in the AFM state, induc- 




tee v) £(eV) 

Figure 3: (left) DOS of a H:ZGNR-4 (both FM and AFM states) in the presence of several periodic 
defects (one defect for 8 ZGNR unit cells): on-top H (top panels), C vacancy (central panels), and 
hydrogenated C vacancy (bottom panels). The local structure of the defects in shown in the middle 
inset. Dashed lines correspond to the pristine ribbon, "negative" DOS is for the minority spin, 
(right) Band structures for both FM and AFM states, black solid lines (red diamonds) correspond 
to the majority (minority) spin. 
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Figure 4: Top panel: structural model of the long ZGNR, used for conductance calculations. Left 
bottom panel: spin-up (G|) and spin-down (Gj) conductance and conductance polarization (Pq) 
of an infinite H:ZGNR-24 in the FM state, whose central section (~197 A) has been randomly 
hydrogenated. Averages (solid thick lines) are taken over 30 random configurations (solid thin 
lines). Shaded (grey) regions represent, at each energy, the standard deviation. Dashed (red) 
curves represent the same quantities for the pristine ribbon. Right bottom panel: same as in (a) for 
the AFM state. 

ing half-semiconductivity.-22»^>2^*^>^ Spin stiffness along the edges, as reported recently,— is 
remarkably high, raising hopes of room-temperature coherence for these one-dimensional systems 
at length-scales comparable to those of microelectronics devices. 

A breakthrough application for graphene ribbons would be in the role of spin valves; indeed, 
energy windows where electronic states with only one spin are available occur under different 



circumstances. We first explore this point in Figure 3 where we compare the DOS and band 
structures of three defected H:ZGNR-4 in the presence of bulk ribbon hydrogenations, carbon va- 
cancies (as e.g. induced by ion irradiation) or both; in all cases the breaking (enhancement) of 
the spin-up/spin-down symmetry (asymmetry) of the AFM (FM) state is observed. On the other 
hand, in any real device no degree of order along the ribbon can be expected, and spin-polarized 
currents along the edges would be overshadowed by the spin-unpolarized transport channels avail- 
able inside the bulk of the ribbon; around Ef, these are dominated by the MLWFs contributions 
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Figure 5: Top panel: spin-up (G|) and spin-down (G|) conductance and conductance polarization 
(Pg) of an infinite H:ZGNR-12 in the FM state, whose central section (~197 A) has been randomly 
hydrogenated. Averages (solid thick lines) are taken over 30 random configurations (solid thin 
lines). Shaded (grey) regions represent, at each energy, the standard deviation. Dashed (red) 
curves represent the same quantities for the pristine ribbon. Bottom panel: same as in (a) for the 
AFM state. 

detailed before. Thus, engineering a ZGNR device for spintronics applications will require remov- 
ing the unpolarized bulk conduction channels while preserving edge conductance. This could be 
achieved with hydrogenations, or ion or electron beams; these are all now recognized as effective 
tools to tailor electronic and transport properties,— >2^^i^>^ and proton irradiation even induces 
ferromagnetism in NM graphitic samples.— '^^'^'^^'^ 

We investigate here this possibility, and show that the availability of spin-polarized edge states 
allows for the design of simple and robust spin-valves. Our strategy is based on the observation 
that most or all functionalizations to the bulk of the ribbon induce sp^ to sp^ re-hybridization, and 
remove the "half-filled" p^ MLWFs from the energy window around £■/?.— This mechanism is 
central to the realization of a spin valve where unpolarized conduction channels are removed from 
the bulk of the ribbon, while preserving at the same time the edge states. Chemical reactions with 



atomic hydrogen represent the simplest route, and we show in Figure 4 the quantum conductance 
calculated for a realistic device: an infinite H:ZGNR-24 (~50 A width) whose central section 
(3840 atoms, ~197 A) has been randomly hydrogenated with a 0.5% defect concentration. The 
conductances G| and G| for both spin channels as well as the conductance polarization Pq (defined 
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as Pg = (G| — Gj) / (G| + Gj)) are averaged over 30 configurations, each with random hydrogena- 



tions. For comparison, we show in Figure 5 the same quantities for the H:ZGNR-12 ribbon (the 
central section now contains 1920 atoms). In the FM state, the spin valve character of such system 
is immediately apparent, with a pronounced spin up or spin down conductance peak just below 
or above Ef, where no significant contribution from the other spin is found. Therefore, tuning 
the applied bias voltage will turn on and off the injection of electrons with a given spin with an 
efficiency that may reach 100% for single devices. In the AFM state, while a random distribution 
of defects would induce a spin polarization due to the breaking of the symmetry between spin up 
and spin down channels, the average effect is zero. Yet, mesoscopic fluctuations mean that for each 
and every device spin valve effects would be very relevant even in this case, as already inferred 
for graphene ribbons with rough edges.— Such effects appear to be more relevant for the smaller 
ribbon ( [Figure 5 ). Irradiation with ion or electron beams, that nowadays can be focused onto areas 



as small as a few nm or a few A in diameter respectively, would lead to similar effects, since the 
creation of vacancies (hydrogenated or not) in the bulk of the ribbon also destroys the unpolarized 



conduction channels (see Figure 3 ). Provided one edge of the ribbon can be physically protected, 
a wide variety of chemical functionalizations - starting from double hydrogenations at the opposite 
edge - would achieve the same objective. 

In conclusion, we have calculated the electronic- structure and quantum conductance in real- 
istic, functionalized graphene nanoribbons with full chemical accuracy, thanks to the use of an 
accurate but minimal MLWFs basis. Defects in spin-polarized ZGNRs - following hydrogenation 
or treatment with ion or electron beams - can be effectively used to remove the unpolarized con- 
duction channels in the bulk of the ribbon, while preserving edge states. Similar results can be 
devised using chemical routes for graphene functionalization/hydrogenation.— For instance, very 
recently the reversible hydrogenation of a graphene sheet has been experimentally proven, lead- 
ing to the sp^ (graphene) to sp^ (graphane) transformation of the carbon network.— The resulting 
asymmetry between the spin-up and spin-down channels makes these systems ideal candidates for 
spin-polarized transport with a very high degree of spin polarization. 
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